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By analyzing trajectories of solid hydrogen tracers in superfluid *He, we identify tens of thou- 
sands of individual reconnection events between quantized vortices. We characterize the dynamics 
by the minimum separation distance 5{t) between the two reconnecting vortices both before and 
after the events. Applying dimensional arguments, this separation has been predicted to behave 
asymptotically as S{t) ~ A {K,\t — to\)^^'^ , where k — h/m is the quantum of circulation. The major 
finding of the experiments and their analysis is strong support for this asymptotic form with n as 
the dominant controlling feature, although there are significant event to event fluctuations. At the 
three-parameter level the dynamics may be about equally well-fit by two modified expressions: (a) 
an arbitrary power-law expression of the form 5(t) — Blt — tol" and (b) a correction-factor expression 
5{t) = A{K,\t — to\)^^'^ (1 + c\t — to\). The measured frequency distribution of a is peaked at the 
predicted value a = 0.5, although the half-height values are a = 0.35 and 0.80 and there is marked 
variation in all fitted quantities. Accepting (b) the amplitude A has mean values of 1.24 ± 0.01 and 
half height values of 0.8 and 1.6 while the c distribution is peaked close to c = with a half-height 
range of —0.9 s~^ to 1.5 s~^. In light of possible physical interpretations we regard the correction- 
factor expression (b), which attributes the observed deviations from the predicted asymptotic form 
to fluctuations in the local environment and in boundary conditions, as best describing our exper- 
imental data. The observed dynamics appear statistically time-reversible, which suggests that an 
effective equilibrium has been established in quantum turbulence on the time scales (< 0.25 s) inves- 
tigated. We discuss the impact of reconnection on velocity statistics in quantum turbulence and, as 
regards classical turbulence, we argue that forms analogous to (b) could well provide an alternative 
interpretation of the observed deviations from Kolmogorov scaling exponents of the longitudinal 
structure functions. 
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I. INTRODUCTION 



Relaxation toward equilibrium requires dissipative pro- 
cesses but can be inhibited by topological defects that 
cannot diffuse. Linear topological defects occur in a va- 
riety of systems such as superfluids p] , liquid crystals 
and superconductors [3J. Dissipation normally accom- 
panies the reconnection of two defect lines that cross, 
change topology by exchanging ends and separate (as il- 
lustrated in Fig. 1). Prime examples T of dissipation 
by reconnection occur in astrophysical plasmas (such as 
solar flares [51[B] and magnetic substorms [315]) and saw- 
tooth crashes in fusion devices [5], where magnetic en- 
ergy is dissipated by the acceleration of nearby particles 
[TCTl im mj. Reconnection has also been studied 

in liquid crystals [2] , superconductors [31 [TSl [T5] , cosmic 
strings [H], viscous [H UHl [2D] and Euler [3T] vortices, 
Bose-Einstein condensates [22] and superfluids [23l [24] 

[23liail3[5Hl[ia[3i[nj[31[MJIMllMl[Ml[S71[3H]. 
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FIG. 1: Schematic of the evolution of reconnecting antiparal- 
lel directed linear topological defects; in the case of vortices 
the arrows indicate the sense of the vorticity. The minimum 
separation between the defects is 5{t) with S{to) = 0. Al- 
though not evident here, reconnection need not be planar 
and, indeed, has been predicted to be intrinsically three- 
dimensional I24ll27l. 



Superfluid "^He behaves as a mixture of two interpen- 
etrating fluids, a viscous normal fluid and an inviscid 
superfluid pLJ. Vorticity in the superfluid component is 
confined to linear topological defects, called quantized 
vortices, which possess circulation values that are inte- 
gral multiples of k — h/m = 9.97 x 10^"' cm^/s, where 
h is Planck's constant and m is the mass of a helium 
atom. Viscosity does not lead to diffuse vorticity in a 
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quantum fluid; this is clearly different than in classical 
fluids. Superfluid "^He driven far from equilibrium con- 
sists of a complex tangle of quantized vortices, known as 
quantum turbulence [23] . Quantized vortex reconnection 
reduces the defect line length and transfers the energy to 
the normal fluid via mutual friction |39l HOI |41] or into 
acoustic and Kelvin wave emission [3TJ [511 ESI 155] , 
This dissipative process is particularly significant in the 
limit of absolute zero temperature owing to the vanishing 
viscosity of the normal fluid component [42l . 

In this article, we characterize experimental data for 
nearly 20,000 quantized vortex reconnection events in su- 
perfluid '*He visualized using micron-sized solid hydrogen 
tracer particles. Previous studies have shown that hydro- 
gen tracers can be trapped by quantized vortices |43l 144] 
in superfluid *He and thereby used to directly visualize 
quantized vortex reconnection [37l [38] . We characterize 
the dynamics of reconnection by measuring the minimum 
separation distance S{t) between the vortices both before 
and after the reconnection event (see Fig. [TJ. Each mea- 
sured separation sequence has been fit to two expressions 
each described by three parameters, namely, an arbitrary 
power-law of the form 

6{t)^B\t~tor, (1) 

and a correction-factor expression 

6{t) = A{K\t-to\)'^^{l + c\t-to\). (2) 

The two expressions suggest distinct physical interpre- 
tations and implications that will be a focus of this 
manuscript. We provide a discussion of relevant previ- 
ous work in Section The details of the experiments 
and our techniques for identifying reconnection events are 
discussed in Section |III| The arbitrary power-law expres- 
sion ([TJ and the correction-factor expression in ([2| are 
discussed and compared in Section |IV[ We present argu- 
ments for statistical time-reversibility and anisotropy in 
Section|V] The effects of reconnection on quantum turbu- 
lence are discussed in Section |VI| An overview of poten- 
tial implications for understanding classical turbulence is 
presented in Section |VII| while Section |VIII| summarizes 
our conclusions and suggests future investigations. 

II. PREVIOUS STUDIES 

Quantized vortex reconnection has been previously 
studied numerically and analytically by employing 
vortex-Hne methods |24l |25l [23 [29||30l [45] and by inte- 
grating the Gross-Pitaevskii equation [51] . To char- 
acterize the evolution of reconnecting vortices, some of 
these theoretical studies examined the minimum separa- 
tion distance d{t) between the vortices as shown in Fig. 
1. Assuming that the only relevant parameter in recon- 
nection dynamics is the quantum of circulation k, dimen- 
sional analysis yields the relation 

6{t)^AiK\t-to\)'^\ (3) 



where to is the reconnection moment for the vortices 
and A is a dimensionless factor of the order unity. One 
may optimistically expect this scaling to be valid for 
length scales between the vortex core size (- 1 A) and 
the typical intervortex spacing (~ 0.1 — 1 mm for our 
work here). Basically, however, ^ should represent 
an asymptotic expression, subject at least to correc- 
tions for longer times. Indeed, slight deviations were 
observed in prior simulations implementing line-vortex 
models [24] [27]. The Gross-Pitaevskii equation is Hamil- 
tonian and time-reversal invariant, but particular solu- 
tions may well break this symmetry. 

Recent experimental studies have demonstrated that 
hydrogen tracer particles may be used to directly visu- 
alize quantized vortex reconnection f37[ I38j. Such tracer 
particles can be trapped on quantized vortices or may 
move with the normal fluid under the influence of Stokes 
drag [15][11[1S]. In [571 [5F] it was found that the form ^ 
must be modified to adequately represent the set of ex- 
perimental data. Specifically, in [57] the separation S{t) 
for t > to for 52 distinct events was fit to the form ([ij 
while in [5S] the alternative form ^ was used both be- 
fore and after approximately 20,000 reconnection events. 
Here we undertake a comparison of these two different 
expressions. 



III. EXPERIMENTS AND TECHNIQUES 
A. Pulsed Counterflow Experiments 

Our experiments are conducted in a cylindrical cryo- 
stat with a 4.5 cm diameter containing liquid '*He. The 
long axis of the channel is vertical with four 1.5 cm win- 
dows separated by 90°. A room-temperature mixture of 
2% H2 and 98% *He by partial-pressure is injected into 
liquid helium [ll] above the superfluid transition tem- 
perature (Tx = 2.17 K). The injected hydrogen freezes 
producing a polydisperse distribution of tracer particles 
with diameters near 1 //m and an initial volume fraction 
of approximately 10^^ — 10~^. We subsequently cool 
the fluid evaporatively to the desired temperature in the 
range 1.70 K < T < 2.05 K. A portion of the hydro- 
gen leaves the observation volume resulting in a volume 
fraction ^ 10~^. The solid hydrogen particles are illumi- 
nated by an argon ion laser sheet that is 8 mm high and 
100 fim thick and has an optical laser power between 2 
and 6 W. A camera gathers 90° scattered light with a 
resolution of 16 fj,m per pixel at either 60, 80, or 100 
frames per second. 

A quantum vortex tangle and its accompanying cas- 
cade of reconnection events is induced by reproducibly 
driving the system away from equilibrium by a thermal 
counterflow [2 [551 I1S]> which increases the total vor- 
tex line length present in the system [25l [39l [40] |4T] . 
Upon cessation of the counterflow, we acquire data 
while the system relaxes toward equilibrium |47j . A 
spiral nichrome wire heater located at the bottom of 
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FIG. 2: Contrast-enhanced negative images of particles 
trapped on reconnecting vortices (top) along with velocity 
vectors (middle) and our interpretation of the pre- and post- 
reconnection configurations of the vortices denoted by the red 
and blue lines (bottom) with time measured from to- The 
green vectors show the background drift that is subtracted 
from all velocity vectors. The red and blue velocity vectors 
(middle) correspond to the red and blue marked vortices in 
the bottom images. The volume fraction of hydrogen in these 
images (10~^) is higher than for all the pulsed counterflow 
experiments discussed below. Reconnection is particularly 
unambiguous in the online movies in [38) . 



the channel 7.5 cm below the observation volume ini- 
tiates the counterflow. A fixed heat flux in the range 
0.064 W/cm^ < q < 0.17 W/cm^ drives the system for 
approximately 5 s, after which it is allowed to relax for 
approximately 10 s before repeating the process. We em- 
ploy a two-dimensional particle-tracking algorithm |48j 
with sub-pixel precision to obtain single particle trajec- 
tories. We estimate the uncertainty in the particle posi- 
tions to be less than 4 /im. While superficially similar, 
this technique is distinct from particle image velocimetry 
(PIV) PIV obtains velocities by computing cross 

correlations of groups of particles, thereby requiring a 
smoothly-varying velocity field. Given the two-fluid na- 
ture of superfluid helium this technique would be unsuit- 
able for the present studies. 



B. Identifying Reconnection Events 

Near the reconnection moment to, reconnecting vor- 
tices move with high, atypical velocities and accelera- 
tions. An example of a reconnection event is shown in 
Fig. [2] (also see Fig. 2 in |37j and note that reconnection 



is clearly evidenced in the online movies in [3S]). In this 
example, the particle density is high so that both vortices 
are marked by multiple trapped hydrogen particles. The 
two vortices merge, exchange tails, then separate as indi- 
cated by the velocity vectors in the middle row of images 
in Fig. I 

Since the hydrogen particles are not completely pas- 
sive, the hydrogen volume fraction in the pulsed counter- 
flow experiments presented here has typically been kept 
one to two orders of magnitude lower than that shown 
in Fig. |2] For such low volume fractions, each identifled 
vortex has only one to a few hydrogen particles trapped, 
thereby minimizing the effects of the hydrogen on the re- 
connection dynamics |46j . A reconnection event is char- 
acterized, then, by a pair of particles rapidly approaching 
or separating. The number of possible particle pairs an- 
alyzed is ~ 10^°, which requires an ad hoc criterion to 
determine likely reconnection events. We deflne particles 
i and j as marking a reconnection event at time t if the 
pairwise separation Sij{t) = \j^i{t) — rj(i)| satisfles 



=,5,,(t±0.25 s)/(5,,(i) >4, 



(4) 



where Vi(t) is the two-dimensional projection of the po- 
sition of particle i at time t and the plus (minus) sign 
indicates particles that separated after (approached be- 
fore) an event, which we label as forward (reverse) events. 
We choose the temporal duration of 0.25 s to allow a suf- 
flcient range to perform the power-law flts to the data 
while curtailing greater times, which are dominated by 
boundary effects and the presence of neighboring vor- 
tices. The criterion ^ excludes all but a fraction of 
possible pairs, namely ~ 5 x 10^ forward and a similar 
number of reverse events. 

It is important to note that we are assuming 6{t) ~ 
Sij{t); however, the particles (z, j) may not be located as 
close as desirable to the point of reconnection. We do not 
observe any correlations between the measured quantities 
discussed below and the initial particle separations or 
the values of as deflned in (|4]); nevertheless, more 
detailed theoretical analyses of vortex reconnection are 
needed to reveal and quantify systematic effects that may 
be caused by interpreting our measurements of dij (t) as 
good approximations to 6{t). 



IV. RECONNECTION DYNAMICS 
A. Arbitrary Power-Law 

We characterize the dynamics of reconnection by mea- 
suring the separation S{t) ~ <^(i) of pairs of particles 
{i, j) that meet the criterion (|4|. As mentioned above, 
previous dimensional and theoretical arguments predict 
that 6(t) behaves asymptotically as a power-law with a 
scaling exponent a = 0.5. To test this hypothesis we fit 
our data to an arbitrary power-law of the form 



(5) 
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FIG. 3: Variation of X^{a) normalized by its minimum value 
Xmin = 0.73 as a function of the scaling exponent a for the 
event shown by the red squares in Fig. |8] below. We choose 
the parameters of the arbitrary power-law fit {a, B, to} that 
minimize 



^ as a function of a defined by i6\ 



The values of B and to are determined by a linear least- 
squares fit of for 500 values of a evenly-spaced 
in the interval < a < 2. For each set of {a, B, to} we 
compute the error in the fit 



1 " 



(6) 



where m denotes the movie frame, tr = 4 /im (0.25 pixels) 
is an estimate of the uncertainty of the particle positions, 
and n = 15, 20, 25 for data collected at 60, 80, or 100 
frames per second, respectively. We then choose the set 
of {a, B, to} that minimizes (see Fig.|3|. 

The measured separations 6(t) for four forward events 
are shown in Fig. [4] along with the predicted asymptotic 
form S{t) ~ {K\t — toiy^"^ for comparison. Fits to ([5| are 
shown as solid lines with the scaling exponent a given in 
the legend. The most frequent fitted exponents cluster 
around the predicted value of a ~ 0.5 and their corre- 
sponding amplitudes B are of order k^^^; however, there 
is a broad spread in both quantities. 

Distributions of a for both the forward and reverse 
events, determined from fifty distinct experimental heat 
pulses, are shown in Fig. p[a). The distributions are 
formed from events with < 4. Approximately 40% 
of the 50,000 pairs that meet the criterion in Q meet 
this criterion. Both distributions are asymmetric but 
peaked within 10% of the predicted value a = 0.5. Fur- 
thermore, as shown in Fig. p[b), events with fitted values 



near 0.5 typically have lower values of x^- 

The amplitudes B for the same events are strongly 
correlated with the scaling exponent a as shown in Fig. 
|6] We find that events with a ~ 0.5 have amplitudes 
B ~ ^/k, as expected from dimensional analysis. How- 
ever, de Waele and Aarts [37] measured B ~ ^ kI^-k in 
numerical simulations of quantized vortex reconnection 
in superfiuid "^He using line-vortex methods; this is ap- 
proximately 30% of our experimentally determined value. 
The time-scales in our experiments differ greatly from 
these numerical simulations; de Waele and Aarts deter- 
mined their value of B for < to — t < 3 /is, whereas our 
time-scales span 1 ms < |t — to| < 100 ms. In addition de 
Waele and Aarts quote an amplitude only for two initially 
antiparallel vortices; other initial orientations might yield 
different values for B. On the other hand, we observe 
only a two-dimensional projection of each reconnection 
event, which would lead us to underestimate i?, poten- 
tially furthering the discrepancy. Clearly, resolving the 
source of this discrepancy warrants additional investiga- 
tion. 

The predicted scaling of a = 0.5 is derived from the 
assumption that the quantum of circulation k is the only 
relevant parameter over the length- and time-scales of 
interest . This assumption is valid in the context of line- 
vortex methods, and holds approximately for the Gross- 
Pitaevskii equation at length scales large compared with 
the core diameter. However, deviations from a = 0.5 
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FIG. 4: Four forward events well fit by the arbitrary power- 
law expression ([5|. Symbols denote the measured separa- 
tion (5(t) of pairs of particles on reconnecting vortices with 
an example error bar cr = 4 /^m while solid lines show fits to 
= B\t — to|° with a given in the legend. The predicted 
asymptotic scaling &(t) = — to|)^^^ is shown by the purple 
dashed line. 



5 





FIG. 5: (a) Normalized frequency distributions of a com- 
puted for 19,150 forward events (black circles) and 18,900 
reverse events (red squares). The mean values of a for for- 
ward and reverse events are 0.68 and 0.69, respectively, (b) 
Two-dimensional contour diagram of versus q for forward 
events. The peak near a — 0.5 with low values of indi- 
cates that ([5| best describes events with dynamics near those 
predicted in ([3|. 



might be obtained, at least conceptually, in two ways. 
First, adapting arguments that were proposed for mul- 
tiscaling solutions of the Euler equation related to in- 
termittency in classical fluid turbulence |50l , one might 
suppose the precise value of k is irrelevant to the dynam- 
ics of reconnection on the length scales we observe; then 
it should be possible to form a continuous family of so- 
lutions with differing values of a. However, it is difficult 
to imagine conditions under which k would be irrelevant 
in our experiments given that the observed velocity mag- 
nitudes seem closely related to k and to the distances 
involved, as expected. 

Second, if another parameter with units different from 
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FIG. 6: Plots of the average fitted amplitude B as correlated 
with the scaling exponent a for 19,150 forward events (black 
circles) and 18,900 reverse events (red squares). The error 
bars indicate the standard deviation of the data within the 
given range of a. The blue diamonds at a = 0.5 show B = 



,1/2 



to compare to the predicted scaling, (a) The entire range 



of q; (b) range restricted to < q < 1. 



K were relevant then it is possible to rationalize recon- 
nection dynamics with a 0.5. If any parameter such as 
a vortex-core length-scale, core surface tension, typical 
intervortex spacing, local velocity gradients, or system 
size, were relevant to the reconnection dynamics, then 
we may construct putative solutions with variable val- 
ues of a. For example, if a core surface tension 7 were 
relevant, we could contemplate an expression of the form 



6{t) = B\t--U 



il/p) 



2q-1 



(7) 



where we use the density p to construct a kinematic sur- 
face tension ^/ p. Note that for a = 0.5 we recapture 
the predicted behavior ([s]). The value of a for a partic- 
ular reconnection event, in this interpretation, is either 
determined by the allowed values of a from the nonlin- 
ear equations of motion, or should that not be unique, 
additionally by the initial and boundary data for each 



particular event. Further deeper theoretical investiga- 
tions are required to determine if such solutions could be 
realized under experimental conditions. 



B. Correction-Factor Expression 

The dynamics of reconnection may alternatively be de- 
scribed by supplementing the predicted asymptotic scal- 
ing of ([3| with a correction factor. The simplest and 
natural expectation is the three-parameter form 



5{t-Q 
(cm) 



5{t)^ A{K\t-t^\)^'\l + c\t~U\). 



(8) 



To test this expression we have performed a linear least- 
squares fit to determine A and c for 500 values of io evenly 
spaced 125 /iS apart (see Fig. I?]) ; we then select the set 
of {A, c, io} that minimizes x defined in 

The measured separations 5{t) for four forward events 
are shown in Fig. [8] along with the predicted scaling of 
^ for comparison. Solid line fits to the correction- factor 
expression in ([s]) describe the data well. Distributions of 
the amplitude A, computed from the same fifty distinct 
experimental heat pulses used to form Figs. |5]and[6j are 
shown in Fig. |9ja). As before, we require < 4. For 
both the forward and reverse events, the distributions 
of A are peaked near unity, in accord with the dimen- 
sional arguments; however as in Fig. |6] the values typ- 
ically exceed unity (so being about three times greater 
than found in the special case studied numerically by de 
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FIG. 7: Variation of X^(^o) normalized by its minimum value 
Xmin = 0.73 (black triangles) and the corresponding correc- 
tion factor c (red circles) as a function of the time origin to 
for the event shown by the red squares in Fig. |8] We choose 
the parameters of the correction-factor expression {A, c, to} 
that minimize defined by ([gJi. 
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FIG. 8: Four forward events well fit by the correction-factor 
expression. Symbols denote the separation 5{t) of pairs of 
particles on reconnecting vortices with an example error bar 
CT = 4 fj.m while solid lines show fits to the correction-factor 
expression ||8|. The predicted asymptotic form S{t) — {K\t — 
to\y^'^ is shown by the purple dashed line. 



Waele and Aarts [21] )• Events with A near unity typi- 
cally have lower values of x^, as shown in Fig.joJ^b), again 
supporting the inferences based on dimensional analysis. 
The distributions of the correction amplitude c for the 



forward and reverse events are shown in Fig. 10 'a). The 
distributions are peaked at c = 0, indicating that many 
events follow rather closely the simple scaling of (|3|. 
However, for both the forward and reverse events, the 
distributions are broad compared to their mean values 
of 0.63 and 0.71 s~^, respectively, signifying strong event 
to event variation. The correction factor also varies sys- 
tematically with the fitted value of the time origin as 
illustrated in Fig. [7] and easily understood algebraically. 
Thus, uncertainties in the estimation of Iq will induce 
corresponding changes in estimates for c. 

It is useful to consider the magnitude of the correc- 
tion term, since larger magnitudes imply a greater de- 
parture from the asymptotic form (|3|. This departure 
can be quantified by |c(t — to)|j where the overbar im- 
plies a time-average over the duration that we use to 
fit the data, namely < |t — io| < 0.25 s. The ensemble 
mean value of \c{t — to) I for both the forward and reverse 
events is {\c{t — to)\) = 0.15, but it ranges from 10~^ to 
values greater than unity. This implies that while the de- 
viations vary from event to event, they typically amount 
to less than ±20 % as evidenced in Fig.[To](b). 

The need to supplement the simple form (|3|, derived 
by dimensional analysis, with a correction factor stems 
from several potential sources. Independent of a spe- 
cific origin, one must always expect a subdominant term 
in asymptotic power-law scaling forms like ([s]) which 
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FIG. 9: (a) Normalized frequency distributions of the ampli- 
tude A for 19,600 forward events (black circles) and 19,300 re- 
verse events (red squares). Both distributions are broad with 
a peak at ^ = 1.25, the means of the forward and reverse 
distributions of A being 1.25 and 1.23, respectively, (b) Two- 
dimensional contour diagram of versus A for the forward 
events. The peak near A — 1 at low values of indicates 
that ([sjl describes optimally events with dynamics close to 
those predicted in Q. 



describe behavior from a micro- or mesoscopic domain 
up to some appropriate infrared cutoff at long times or 
large length-scales. A functional form including a cor- 
rection factor along with the dominant power-law allows 
for a crossover between scales. In our case, one certainly 
should expect deviations from ^ on length scales compa- 
rable to the typical intervortex spacing of 0.1 to 1.0 mm. 
Indeed, if one introduces a correction length scale, say 
by rewriting the correction term as 

c\t-to\ = ±K\t-to\/f, (9) 

for forward and reverse trajectories, one finds that the 
corresponding forward and reverse mean values of c cor- 
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FIG. 10: (a) Normalized frequency distributions of the cor- 
rection amplitude c for 19,600 forward events (black circles) 
and 19,300 reverse events (red triangles), (b) Normalized dis- 
tributions of the magnitude of the correction factor \c{t — to)l 
time-averaged over 0.25 s for each trajectory for the forward 
(black circles) and reverse (red triangles) events. 



respond to I — 0.40 mm and 0.38 mm, respectively. Thus 
the dominant correction may well represent the influence 
of neighboring vortices and their ability to distort the 
observed trajectories. 

In addition, however, other spatial and temporal as- 
pects of the local environment may also significantly af- 
fect the dynamics of reconnection, beyond the leading 
behavior accounted for by dimensional analysis. Local 
velocity gradients or other initial and boundary condi- 
tions could all produce deviations from pure square-root 
scaling, thereby necessitating a nonzero correction factor. 
More intriguing theoretically, however, and challenging 
experimentally, is the possibility of nonanalytic correc- 
tion terms such as cg\t — ioT with 9 nonintegral. The 
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presence of such terms with nontrivial values of ~ 0.5 
is well established in the study of critical phenomena in 
ferromagnets, superfluids, at gas-liquid transitions, etc.: 
see, e.g., [ST1|511[53]. 



C. Comparative Assessment 

The two expressions ^ and ^ are both modified ver- 
sions of the asymptotic dynamics which suggest some- 
what distinct theoretical interpretations. In both cases, 
the majority of the fitted events exhibit behavior very 
similar to the dimensional predictions (i.e., a = 0.5 or 
c ~ 0); although, both also show strong variations from 
event to event. One might hope to distinguish the quality 
of the two fits by comparing the observed distributions 
of x^i the overall deviation in the fits. In fact and un- 
surprisingly, the two distributions shown in Fig. [TT] are 
rather similar and no firm basis for making any distinc- 
tions emerges. 

On balance at this point we favor the correction-factor 
expression as best embodying our experimental data for 
quantized vortex reconnection. Our typical observations 
of the dynamics show only relatively slight deviations 
from those predicted by dimensional analysis. It there- 
fore appears that the dominant parameter is indeed the 
quantum of circulation k, which sets the leading scaling 
exponent of a = 0.5. We expect the deviations from 
the corresponding asymptotic form in our experimental 
range to be caused by the local environment and the ini- 
tial and boundary conditions of the event, as opposed to 
other parameters such as a surface tension of the vortex 




FIG. 11: Normalized frequency distributions of from both 
the forward and reverse events fit to the arbitrary power-law 
expression (black circles) and the correction-factor expression 
(red triangles). 



cores. Indeed, such parameters would likely vary with 
temperature and we have not observed any correlations 
between our fit parameters and the temperature of the 
system over the range 1.70 K< T < 2.05 K. Experiments 
and numerical simulations that control the local environ- 
ment (velocity gradients, neighboring vortices, strains, 
etc.) as well as the initial and boundary conditions (con- 
figuration of the vortices, initial velocities and curvatures, 
etc.) could directly test this hypothesis and are clearly 
desirable. 



V. TIME-REVERSIBILITY AND ANISOTROPY 

The Gross-Pit aevskii equation as used for quantized 
vortex reconnection [26] 1341 is fully symmetric under 
time reversal; thus solutions of the equation may also 
be time-reversible symmetric. However, many previous 
theoretical works have concluded that reconnection dis- 
sipates energy by emitting acoustic and Kelvin waves 
[3TJ [32l |33l |35] that may be absorbed by the boundaries, 
and so would break locally the time-reversibility. We 
have compared the pre- and post-reconnection dynamics 
by separately fitting forward and reverse events. All of 
the distributions of the fit parameters (see Figs, [s] 6] [o] 



and 



10 1 for the forward and reverse events show stri 



kmg 



similarity, as would be expected if the dynamics were sta- 
tistically time-reversible. While it is clear, however, that 
some energy is dissipated overall in our experiments, as 
evidenced by the decay of the turbulent state, it is not 
evidenced in the statistics of individual events. 

In addition to the fit parameters, we may also inves- 
tigate the total displacement of the vortices before and 
after events. We define the displacement vector of parti- 
cle i as 

Avi = ri{to + 0.25 s) - ri{ta) = xAxi + zAz^. (10) 



Figs. 12 'a) and (b) show the x- and z-components of 
the displacement vectors for all of the particles identi- 
fied with forward (black) and reverse (red) reconnection 
events. Indeed, the forward and reverse displacement 
vectors also show striking similarities. The displacement 
vectors appear weakly correlated in the a;-direction and 
anti-correlated in the z-direction. This anisotropy is 



clearly exhibited in Fig. 12 'c), which shows the differ- 
ence in the displacement vectors of the pairs of particles 
on reconnecting vortices. The forward and reverse events 
are found to be equally affected by the these anisotropic 
effects. 

We believe the anisotropy arises from the polarizing 
effect of the z-directed, initiating thermal counterflow. 
A possible interpretation of the anti-correlation in the z- 
direction is that the vortices are typically aligned or anti- 
aligned with the direction of the counterflow. Previous 
studies have also observed |54] or argued for [55l |56l |57] 
the presence of anisotropy in counterflow turbulence; al- 
though, subject to interpretation, it is not clear that we 
agree on the polarization of the anisotropy. 
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FIG. 12: Scatter plots of (a) the i-component and (b) the ^-component of the displacement vector Ar^ = ri{to +0.25 s)— ri(io) 
for all forward (black) and reverse (red) event pairs (labeled as i = 1 and j = 2 for all pairs), (c) Normalized frequency 
distributions of the differences in i-components (black circles) and the i-components (red triangles) of the displacement vectors 
for each pair of particles on reconnecting vortices. The central peak at Axi — Aa;2 = signifies a strong correlation of the 
particle trajectories in the a;-direction while the minimum at Azi — Az-^ — implies a strong anti-correlation in the ^-direction 
associated with the z-directed thermal counterflow. 



Overall the close statistical similarity of the forward 
and reverse events suggests an effective equilibrium has 
been established in quantum turbulence on the time 
scales (< 0.25 s) we have investigated. We have only 
limited statistics to directly compare the dynamics be- 
fore and after a given event, the details of which will be 
reserved for a future publication. However, they are suffi- 
cient to show that individual events are not time-reversal 
invariant. An important future direction for both exper- 
iments and numerical simulations entails understanding 
the interplay between irreversible individual events and 
the reversible statistics of quantized vortex reconnection. 



VI. EFFECTS OF RECONNECTION ON 
QUANTUM TURBULENCE 

Reconnection has long been considered to play an im- 
portant dissipative role in quantum turbulence. Vinen 
[39t l40l |4T] described how the balance of reconnection 
and mutual friction leads to saturated vortex line lengths 
in counterflow turbulence, which is analogous to the sat- 
uration of dynamo action [58 or the magnetorotational 
instability |i59] produced by magnetic reconnection in as- 
trophysical plasmas. A great deal of recent research, 
though, has focussed on the behavior of quantum tur- 
bulence on length scales sufficiently large that the in- 
teractions of individual quantized vortices may be ne- 
glected. These previous works concluded that on such 
length scales quantum turbulence shares many charac- 
teristics with classical turbulence [17, 60, Gil IS^ 1551 [Ml 
IMllMlEIllMllMlllQllITllIllllSllTitJSi^ However, 
the strength of the evidence has been questioned |77j . 

The assumptions used to argue for the classical na- 
ture of quantum turbulence break down on length scales 



smaller than the typical intervortex spacing, which is 
what we have probed with our measurements. It is 
clear from the online movies (see [38]) that our turbu- 
lent states in superfluid ''He differ drastically from those 
observed in classical fluids as a result of the topological 
interactions of the quantized vortices. Speciflcally, recon- 
nection produces anomalously large velocities in highly- 
localized areas, which are not diffusively smoothed by 
viscosity. If length-scales are evolving asymptotically as 
R{t) — A\K{t — to)r^^, then we expect the velocities to 
scale as 



vit) 



(11) 



which far exceed typical fluid velocities when t — > to- 
Note, however, that we expect such velocities to be cut- 
off by the speed of first sound. 

To model the velocity statistics of a turbulent state 
in superfluid "^He characterized by many reconnection 
events we may use the transformation 



Fr,{v) - \dt/dv\Frt[t{v)], 



(12) 



where T'Ty{v)dv is the probability of observing a velocity 
between v and v + dv at any time while T'Tt{t)dt is the 
uniform probability of taking a measurement at a time 
between t and t + dt. Accepting the relation (111, we 
predict for large v (small |i — to|) the behavior |38j 



Pr„(?;) cx \dt / dv\ cx |t;| 



(13) 



A probability distribution function of the velocities de- 
rived from all tracer trajectories (not only those identi- 
fied as marking reconnection) is shown in Fig. |13| The 



velocities are computed by forward differences of the par- 
ticle trajectories. A single-parameter fit of the form 
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FIG. 13: Frequency distribution function of tracer particle ve- 
locities resulting from the decay of a turbulent state produced 
by a 0.17 W/cm^ heat flux at a temperature T = 1.90 K. The 
distribution is scaled by the standard deviation ~ 0.08 
cm/s. The straight line (in red) is a fit to Pr„(ii) = an""^ for 
v/ct„ > 1.5. 



Pr^(t;) — av where v 
solid line for comparison 



— {vl + v'^y^^, is shown as a 



Evidently, the simple argu- 
ments used to derive (13) are able to predict the tails 



of the velocity distributions in reconnection-dominated 
quantum turbulence. 

The velocity statistics of classical turbulence 78J in 
both experiments [75] and numerical simulations [80, 81 
are nearly Gaussian over several orders of magnitude in 
probability. Such statistics are in stark contrast to the 
power-law statistics found in Fig. [TSj We attribute this 
distinction to the topological interactions of the quan- 
tized vortices, which do not exist in classical turbulence 
where the velocity field is diffusively smoothed by viscos- 
ity. It is important to note that the normal fluid present 
in our experiments is relatively quiescent and potential 
future directions could include examining the velocity 
statistics for the case where both the normal fluid and su- 
perfluid are turbulent, since the two fluids couple through 
friction acting on the quantized vortices [33 HOI IH] ■ 



VII. IMPLICATIONS FOR CLASSICAL 
TURBULENCE 

The scaling properties of velocity correlations in clas- 
sical turbulence have received a great deal of attention. 
In particular, much debate has addressed the values of 
the exponents C„ of the longitudinal structure functions 



r''" , where Aur = w(x + r) — u{x) for a single 



In analogy to the predictions of a = 0.5 for quantized 
vortex reconnection, Kolmogorov used dimensional ar- 
guments to predict = «/3 [SI]. His theory pivots on 
the assumption that the dissipation per unit mass e is 
the only relevant parameter in the observed correlations 
and spectra. However, experimental observations report 
values of C„ that deviate slightly from Kolmogorov scal- 
ing [85] IHIl [Ml IMl IM] ■ Typically, arbitrary values of the 
exponents Cn are fit to the data [7^1 |5U l |55 1 |Hi l |55 1 |M l 
|871[88llig[9ni[Ml[9aiMl[Ml[9a[Ml[^, in analogy to 
our arbitrary power-law expression in ([Sl. It should be 
noted, however, that models with variable C„ presuppose 
that arbitrary exponents are allowed, based either on the 
irrelevance of e, or on dimensional grounds, by the ad- 
mission of other relevant quantities that yield new power 
law forms [as illustrated in ([t])]. 

We argue, though, that another option is available - 
that of correction-factors representing subdominant scal- 
ings. Similar connections between critical phenomena 
and turbulence have been explored previously [SH llOOj . 
Here we specifically suggest that individual events in 
classical turbulence might be modeled both by a dom- 
inant Kolmogorov term and a correction factor arising 
from various causes including at least the local neighbor- 
hood conditions and finite-size effects. The basic Kol- 
mogorov scaling derives from the fact that with vanish- 



ing viscosity one obtains (Au^) 



ler. If one inter- 



velocity component u parallel to r [7S] [Bni 82, 83, 84 

[S3iMi[E3iEi[Ei[nniinii[ni[niiMi[n3[Mii9L^,99j. 



prets these statistics as stemming from many individual 
"Kolmogorov events" with Aur ~ S/{t — t^), then sub- 
stituting r — > (5 above yields, at least on a dimensional 
basis, S ^ e^/^(t — to)"^^^, which is also known as Richard- 
son scaling jlOlj . This can be obtained directly from a 
dimensional argument if S depends only on e and time. 
Note the units of e are rn^/s^. If one then extends this 
model to a correction- factor expression, similar to (|8|, 
with an appropriately chosen correction, one might ob- 
tain behavior that would be difficult to distinguish from 
fluctuating power laws, though with a rather different 
interpretation. 



VIII. CONCLUSIONS 

In conclusion, we have observed the dynamics of in- 
dividual reconnection events in superfluid "^He and their 
effects on decaying quantum turbulence. Although we 
observe significant deviations that vary from event to 
event from the mean behavior, the typical dynamics are 
close to those predicted by dimensional arguments. We 
regard this as our major finding. The deviations may be 
accounted for in two separate ways: (a) by supposing the 
scaling exponent of the dynamics can fluctuate as in (|5|, 
or (b) by recognizing that the predicted power laws must 
be supplemented by a correction factor such as in (|8|. 
The two three-parameter expressions describe the data 
almost equally well from a perspective, but suggest 
distinct physical interpretations. Thus, we observe that 
a variable scaling exponent should result from either a 
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lack of importance of the precise value of the quantum 
of circulation k or from the competing relevance of an- 
other physical quantity of distinct dimensions (such as a 
length-scale, surface tension, etc.). 

On the other hand and more naturally, we interpret 
the correction factor as arising from initial conditions 
and boundary effects, such as the vorticity distribution 
and intervortex spacing, and from properties of the local 
environment at reconnection, such as velocity gradients, 
pressure gradients, and thermal fluctuations. Since the 
dynamics appear to be well characterized by the pre- 
dictions that assume that the only relevant physical pa- 
rameter is the quantum of circulation k, we believe our 
data indicate that the environment, as opposed to other 
parameters, is most likely the origin of the observed de- 



viations. Further investigations experimentally and the- 
oretically could focus on; (i) determining if each recon- 
nection event is time-reversal symmetric, {ii) considering 
alternate forms of the correction term, such as co\t — to\^ 
with 6^1, and {in) systematically changing the initial 
and boundary conditions as well as the local environment 
near reconnection and investigating the resulting devia- 
tions from the dimensionally predicted asymptotic form. 

We thank Makoto Tsubota, Carlo Barenghi, Joseph 
Vinen, Nigel Goldenfeld, Christopher Lobb, Marc Swis- 
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terials at the University of Maryland. 
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